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Abstract 

It is shown how nucleon-nucleon potentials can be defined in N dimensions, 
using dimensional regularization to continue amplitudes. This provides an easy 
way to separate out contact ((5-function) terms arising from renormalization. 
An example is worked out several ways for the case of two scalar particles 
exchanged between nucleons, which involves a very simple loop calculation. 
This leads to a Feynman-parameterized representation for the nucleon-nucleon 
potential. Alternately, a dispersion representation can be developed leading to 
a different, though equivalent, form. 



It might be surmised that there exist no new ways to develop nuclear potentials 
after decades of experimentation with techniques. We present below a variation on 
several such techniques, which, to the best of our knowledge, has not been used before. 
As an illustration, we perform in several different ways a simple calculation involving 
a closed meson loop, which leads to a two-boson-exchange potential (TBEP). 



Figure 1: The nucleon-nucleon potential resulting from the simultaneous exchange of 
two scalar particles, "a" and "b" . 

Figure 1 shows the interaction of two nucleons "1" and "2" via the exchange of 
two different scalar mesons, "a" and "b" . The corresponding Lagrangian is 

LjNT = >^NNab. (1) 
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The nonrelativistic potential corresponding to Figure (1) can be developed in old- 
fashioned second-order perturbation theory in the usual way[|l| for a nucleon-nucleon 
separation, r: 

''^""^ J (277)3 ^2Ea) J (27r)3 {2Eb){Ea + Et,)' ^ ' 

This awkward and inelegant expression is nevertheless very easy to interpret. The 
initial factor of two accounts for two ways the mesons can be exchanged: "1" to "2" 
or from "2" to "1" . The factors {2Ea) and {2Eb) are the wave function normalization 
factors for "a" and "b" {Ex = {ql+ml)^), while —{Ea+Eb) is the energy denominator 
(neglecting nucleon recoil) and the integration is over the mesons' phase spaces. This 
can be cast into a more conventional form by using0 

EaEb{Ea + Eb) nJo {Ej + (3^){Ei + P^) ' 

which allows a convolution representation to be written in terms of separate Yukawa 
functions for a and b exchange: 



2A^ f°° dye {y'^ + ml - ml){y'^ + ml - ml) 



(46) 



r - 2y\ml + ml) + {ml - ml)']^ . {Ac) 

Jnia+mh V 



(47r)3r2 Jma+mt, y^ [y4 _ 2^/2(^2 _)_ _^ (^2 _ ^2j 

-2A^ r°° dy _^ 
(47r)3r Jma+mt y 

We have made a change of variables to obtain eq. (4b) and an integration by parts to 
realize eq. (4c). Our points can be made if we further assume that one particle mass 
vanishes or that = m^ = m; the latter leads to an elementary integral which we 
evaluate: 

y 



Vu{r) = -—^ dye-y^ " (5a) 

{4:71 j'^r'' J2m (y2 _ 4^2^ 2 

— 4mA^ 

Note the factor of (47r)^, one power from the configuration space propagator and two 
powers from the loop integral. Both are necessary for dimensional power counting to 
work0. 

On the other hand, we can calculate the Feynman diagram associated with the loop 
in Figure (1). It is easy to show|Q that the potential equivalent to that amplitude, 
M, should have an additional factor of z: V = iM . Evaluating the diagram we find 

M{q) = {N,N,){N2N2)X'l{q') , (6a) 
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where q is the transferred momentum. The divergent loop integral is given by 

2 f d^k 1 
^ " y (27r)4 (P - m2)((g - A;)^ - ml) ' ^^^^ 

and dropping the unimportant nucleon spinor factors we find 

Viq') = t\'l{q') . (6c) 

Evaluating divergent loop integrals requires regularization, and for a variety of rea- 
sons, dimensional regularization is the method of choice H today. This yields 



(27r)^(A;2-m2)((g-fc)2-m2^ ' 



where the number of dimensions has been extended from 4 to N, and for N < 4 the 
integral is convergent. The renormalization scale /i keeps the overall dimensionality 
of the integral the same. 



Using the Feynman parameterization 

1 _ /•! dz 
ah Jo [a + (b — a)z]'^ ' 

and shifting the variable to k' = k — qz leads to 



(7) 



/jv(g ) = dz j ^2vr)^ - [ml - [ml - ml)z - q^z{l - z))]^ ^^""^ 
= ^(47r/i2)2-f r(2 - dz [ml - {ml - ml)z - q^z{l - z)]f . (gj,) 

The latter result is a standard form for one-loop amplitudes^. Usually at this point 
one writes 4 — iV = e and performs an e-expansion, leading to a divergent constant 
(which generates a contact term, (NN)"^, in the nucleon-nucleon force) and a finite 
logarithm. We eschew this approach and define a potential in space-time or n 
space dimensions {n = N — 1): 

- / ^e'^n^A^/^(a] , (9) 

where we choose to work in the frame where qo = {In is a function of q^ = q^ — q^)- 
Alternatively, the energy transfer, go? is small and of order {^Y (i.e., a relativistic 
correction) and can be dropped. Both q and r are vectors in n space dimensions. 
The angular integrals can be evaluated which leads to 

f "2" [- -1 

VN{r) = -—ir dqqtjn,(qr) tX^I^i-q^) , (10) 



(27r)t Jo 
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which can be verified easily for n = 3. Note that we are using = as an inte- 
gration variable, which accounts for the sign change in the argument of In- Inserting 
expression (8b) for 1^ and performing the g-integral leads to: 



^^^'^ " v^(4.)%ir"-f io [.(l-.)]-f ' ^''"^ 

P\z) = ^ > . (116) 

z[i — z) 

The factors of 2 and vr clearly depend on n (or A^). Note that the divergent (for 
= 4) F-function has disappeared and the result is finite for r ^ 0. Standard Bessel 
function identities]^] allow us to write for = 4: 

— 2A2 /•! 

= — — - / dz e-^"(l + f3r) . (12) 



where 



(^4-7]-) 37-3 Jq 

For rria = rrih this can be shown |Q] to equal eq. (5b). Note that this representation 
for the potential involves the Feynman parameterization variable, z. Our derivation 
(12) is both elegant and more directly related to the Feynman loop diagrams than 
is the conventional derivation. Indeed, dimensional regularization was developed for 
these diagrams p 

The reason why the divergent factor r(2 — y) disappears is that we implicitly 
renormalized the loop graph when we developed the potential. Performing an e- 
expansion (r(2 — y) ~ f + finite) in eq. (8b) leads to a contact term (~ S^{r)). By 
keeping the internucleon separation, r, finite, these terms don't arise. They have been 
"regularized" away. Note that V^lr) 

^ 2jv— 5 for small v and is progressively more 

singular for increasing A^. 

We can derive this result and an alternative representation (equivalent to eq. (4c), 
for A^ = 4) by using the analytic properties of the field-theoretic amplitude. In the 
t-channel of Fig. (1) the reaction N + N a + 6is open, and this information 
allows construction of a potential, an old but very useful technique[|]. The amplitude 
JAr(g^) in eq. (9) is real for < 0, but for > it develops an imaginary part 
above the threshold for the reaction, which requires > {ma + mh). Switching to the 
integration variable = used in eq. (10), this implies branch cuts for < or 
along the imaginary axis in the complex g-plane indicated in Fig. (2). The left-hand 
cut arises from g§J^_i(ga;). Writing Ju{z) = | H^\z) + HI)^\z) and noting that 
Hj)^^ behaves asymptotically as e*^ while Hj)^^ behaves as e~*^, we can write 

Jn = dqq'^J^_,{qr)lN{-q') = \ J^dq q^ Hf_,{qr)lN{-q') , (13) 
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Figure 2: The analytic q-plane illustrating branch cuts discussed in the text. The 
integration contour-C, which defines the Fourier transform, can be deformed to C, 
which gives the dispersion representation for the potential. 

usmg ^a)(e'^r) = i/(2)(^^-)g-7ri(i.+i)_ rpj^g integral along C can be continuously de- 
formed into C because the integral vanishes exponentially in the upper half of the q- 
plane. Using the properties of the H^P function for imaginary argument and q — ix±e, 
we then find the simple and elegant result: 



which provides an excellent representation for the potential. Note that the range of 
the force varies from {nia + mi,) to oo. Putting everything together from eqns. (8b) 
and (10), we obtain 



where (3'^ — [ml^ — (m^ — ml)z]/z{l — z). Evaluating the imaginary part leads to 




(14) 





(15a) 





N 1 



(156) 
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where the singularity at = 4 has disappeared. 



Two options are available for further simplification. One can first evaluate the x 
integral, which leads immediately to eq. (11a). Alternatively, one can perform the 
z integral first. The function i3{z) is larger than x^ for values of z near and 1, 
making the ^-function vanish. The argument of that function is positive between two 

values of z [z± = ^^^^ ± \ [(1+ -^f), where Am^ = - ml which 

resets the limits on the integral. The resulting integral is a beta-function [7], and one 
obtains 



This reduces to eq. (4c) for N — A. 



(16) 



Equations (11) and (16) are our principal results and are equivalent. The former 
is a potential defined in N dimensions which exploits the Feynman parameterization 
of the loop integral. The latter is a dispersion representation which exploits the 
analytic properties of the amplitude, and was also derived using a trick and simple, 
old-fashioned perturbation theory (in 4 dimensions). Infinities which naturally arise 
are regularized and disappear from the final form. Dimensional regularization was 
used to define the continuation to N dimensions, and both our approach and our 
techniques maintain a strong connection between the potential and the underlying 
field theory. 
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